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The pion form factor is calculated in the framework of the charge-conjngation invariant Covariant 
Spectator Theory. This formalism is established in Minkowski space and the calcnlation is set up 
in momentum space. In a previous calculation we included only the leading pole coming from the 
spectator qnark (referred to as the relativistic impnlse approximation). In this paper we also include 
the contributions from the poles of the quark which interacts with the photon and average over all 
poles in both the upper and lower half planes in order to preserve charge conjugation invariance 
(referred to as the C-symmetric complete impulse approximation). We find that for small pion mass 
these contributions are significant at all values of the four-momentum transfer but, surprisingly, 
do not alter the shape obtained from the spectator poles alone. 
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I. INTRODUCTION 

With the 12 GeV accelerator upgrade at Jefferson Lab, 
the charged pion form factor EV will be known with high 
precision up to momentum transfer « 6 GeV^ [T]. 
This measurement will cover the interesting region where 
the product as a function of Q^, reaches a maxi¬ 

mum and afterwards flattens down, and will help resolve 
the current discrepancy between the results for the 717*7 
transition form factor obtained by the Babar and Belle 
Collaborations. 

In combination with theoretical calculations, the new 
measurements will narrow the uncertainty about the 
smallest momentum transfer at which the description 
based on asymptotic parton distribution functions is still 
valid. One may also certainly expect that the forth¬ 
coming data will clarify the mismatch between physical 
reality and perturbative QCD predictions in the region 
around ~ 6 GeV^ m- This paper focuses on the pion 
form factor in the small region that is planned to be 
covered by the new experiment. Part of the interest on 
this region lies also in its vicinity to the timelike sector 
and on the extrapolation of the derivative of the form 
factor to < 0. The knowledge of the pion form factor 
enters into the evaluation of baryon form factors in the 
region near « 0, and its behavior in the timelike re¬ 
gion helps the interpretation of dilepton production data 
from heavy ion collisions. 

Theoretically, the pion is the consequence both of the 
non-perturbative character of a quark-antiquark bound 
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state, and of Spontaneous Chiral-Symmetry Breaking 
(SySB). Its role in nuclear structure and nuclear dynam¬ 
ics is crucial. The pion cloud is seen to contribute to the 
structure of the nucleon and its excitations, through the 
coupling to external photons. Also, the exchange of pions 
between nucleons dominates their interaction at large dis¬ 
tances and is the primary origin of the tensor force that 
has a decisive influence on the structure of nuclei. 

The non-perturbative dynamics of the pion and other 
hadronic systems has been addressed by constituent 
quark models [SE] and QCD sum rules. These ap¬ 
proaches do not provide a comprehensive and global 
view for both light and heavy mesons, and baryons, and 
they cannot avoid a delicate fine-tuning between many 
parameters. More recently, QCD simulations on the 
lattice mM, light-front formulations of quantum field 
theory EHn, as well as models based on the Dyson- 
Schwinger approach and mass gap equation [I2HI^ , start 
to provide an integrated account of mesons and baryons. 

In this paper we use the Covariant Spectator The¬ 
ory (GST). In common with the Dyson-Schwinger frame¬ 
work, it generates a dynamical quark mass, which is 
a function of the momentum, and this dressed mass is 
consistent with the two-body quark-antiquark dynam¬ 
ics. Although the GST equations share with lattice QCD 
and Dyson-Schwinger equations this important dynam¬ 
ical consistency, in contrast to those approaches GST 
equations are solved in Minkowski space. Therefore in 
this formalism the extension of results from the spacelike 
to the timelike region does not imply further work or the 
use of a different representation. 

In Ref. [50] we presented the first calculation of the 
pion form factor based on the GST Bethe-Salpeter equa¬ 
tion (CST-BSE) and the GST Dyson equation (CST- 
DE), using a dressed quark mass function calibrated to 
fit existing lattice QCD data [H]. In the present pa- 
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per, as well as in Ref. [20], the CST interaction kernel 
in momentum space has the form of a (5-function plus a 
covariant generalization of the linear confining interac¬ 
tion. The fact that this model satisfies chiral symme¬ 
try is worth emphasizing; this is ensured by choosing a 
relativistic generalization of the confining interaction Vl 
which decouples, in the chiral limit, from both the one- 
body CST-DE and from the two-body CST-BSE. In the 
calculation of the pion form factor in Ref. (201 we used the 
relativistic impulse approximation (RIA). However, we 
know this approximation will break down for small pion 
masses, since under these conditions the pole contribu¬ 
tions of the struck quark, which are omitted in RIA, are 
not negligible [2Q]. The contribution of these additional 
poles to the charged pion form factor are calculated here, 
completing the results obtained from the RIA. In addi¬ 
tion, we average over all of the propagator poles in both 
the upper and lower half planes, making our calculation 
consistent with the charge-conjugation symmetric equa¬ 
tions from which the pion vertex must be calculated m- 
This improvement is referred to as the C-symmetric om- 
plete impulse approximation (C-CIA). A still more exact 
result can be obtained by adding interaction currents, 
a dynamically calculated pion vertex function, and the 
full dressing of the quark current to the C-CIA. These 
are planned for future work (for more discussion, see the 
final section). 

This paper is organized as follows: In Section [T^ a brief 
review of the CST formalism is given. In Section |III| we 
present the ingredients for the calculation of the triangle 
diagram, and in Section jlV] the formulas for the contribu¬ 
tions of all the poles to the pion form factor. In SectionfVj 
we present the results and compare to the ones obtained 
with RIA. Finally, Section [VT| presents a short summary 
and conclusions. 


II. BRIEF REVIEW OF THE ONE- AND 
TWO-BODY CST EQUATIONS 

In this section we briefly describe how the CST is ap¬ 
plied to the description of quark-antiquark mesons. For 


technical details, and references to earlier work, see I2D1- 

ESj. 

In the four-dimensional BSE [25] for heavy-light 
mesons, it is known m that cancellations occur between 
iterations of ladder diagrams and higher-order crossed- 
ladder diagrams in the complete kernel. The omission of 
crossed-ladder diagrams and part of the pole contribu¬ 
tions of the ladder diagrams from the kernel can there¬ 
fore give a better approximation to the exact BSE than 
the ladder approximation does. Also, efficiency in the 
summation of the series can be gained by recognizing 
these partial cancellations already at the outset. The 
fundamental idea of CST is then to reorganize the Bethe- 
Salpeter series into an equivalent form—the CST equa¬ 
tion. This leads to a redefinition of both the complete ker¬ 
nel and the (off-mass-shell) two-particle propagators in 
the intermediate states. The resulting three-dimensional 
equation, the one-channel CST (or Gross) Bethe-Salpeter 
equation [2B], CST-BSE for short, is manifestly covari¬ 
ant. An important feature is that—unlike the BSE in 
ladder approximation—the CST-BSE has a smooth non- 
relativistic limit, defining a natural covariant extension 
of the quantum mechanical Dirac and Schrodinger equa¬ 
tions to quantum held theory. 

In the heavy-light mass case, the redehned propaga¬ 
tors are obtained by keeping only the positive-energy pole 
contribution from the heavy quark propagator in the en¬ 
ergy loop integration. The heavy quark is then on its 
positive-energy mass shell. 

In the case of two light quarks, the CST-BSE includes 
an explicit charge-conjugation symmetrization. The ver¬ 
tex functions of 7 r+ and tt~ are connected by charge 
conjugation and, therefore, both positive- and negative- 
energy quark poles must be included. This is realized by 
effectively averaging over the sums of all quark poles in 
the lower and in the upper complex-energy half plane, 
which generates the charge-conjugation-symmetric CST- 

BSE [211 [23, 


r(pi,P 2 ) = -l:Zo f [v(p, k - P/2)A(fc)r(fc, k - P)Sik -P) + V(p, k + P/2)S{k + P)T{k + P, k)K{k) 

^ Jk L 

+V{p, -k - P/2)A{-k)T{-k, -k - P)S{-k - P) + V{p, -k + P/2)S{-k + P)r{-k + P, -fc)A(-fc) 

= i / V{p,k)S{k + P/2)r{K + P/2, K- P/2) S{k- P/2), (1) 

JkO 


where the three-dimensional covariant integration vol¬ 
ume element is 


d^k m 

Wf^k' 


( 2 ) 


and 



ko propagator 
poles only 


- E 

2 

propagator 
pole terms 



( 3 ) 
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In Eq. Q T{pi,P 2 ) is the (4 x 4) bound-state vertex 
function with pi = p ^ and —p 2 = —p + ^ the four- 
momenta of the outgoing quark and antiquark (respec¬ 
tively); ki = K + ^ and —k 2 = —k + ^ are the interme¬ 
diate four-momenta for the quark and antiquark (respec¬ 
tively); P is the total bound-state momentum, p and k 
are the relative momenta; k = {Ek, k) is the on-shell four- 
momentum with Efc = y/m? + k2; V{p, k) = V{p,k;P) is 
the interaction kernel with the general structure 

V(p,K)A’ = ^V-(p,K)0,dfO,, (4) 


and Zq = Z{w?). For E(p) = 0, S(p) becomes the bare 
propagator denoted as So{p). For the models we are us¬ 
ing, Bip^) = 0 and Z(p2) = l. 



where the sum i = {S', P, V,A,T} is over the five possible 
invariant structures that could contribute: scalar, pseu¬ 
doscalar, vector, axial-vector, and tensor. The dressed 
quark propagator, S(fc), projection operator, A(fc), and 
the numerator function, iV(fc), are given by 


S{k) 

m 


1 

mo — -f Z,{k) — ie ’ 
M(fc2) + 11: _ N{k) 
2M(fc2) “ 2M(fc2) ’ 


(5) 

( 6 ) 


where M{k‘^) the dressed quark mass function, mg the 
bare quark mass, and E(fc) is the quark self-energy. For 
consistency, F(fc) is the solution of the one-body CST- 
DE using the same interaction kernel V that dresses the 
quark-antiquark vertex. One obtains [21] 


S(p) = ^ {V{p, fc)A(fc) + V(p, -fc)A(-fc)} 

= -i[ V{p,k)S{k). (7) 

JkO 

The self-energy is written 

T.{p) = A{p^)+fB{p^) (8) 


and then 


S{p) = Z{p^ 




M^{p'^) — _ ig 


= Z{p^) 


D{p^y 


(9) 


where D{p^) = —p^ — ie is the denominator of the 

propagator, and the mass function M(p^) and the wave 
function normalization Z(p^) are 


M(p2) 


A{py + Too 
1-P(p2) ’ 
1 

1 - P(p2) ’ 




FIG. 1. (Color online) The two triangle diagrams for the 
electromagnetic pion form factor. The top diagram describes 
the interaction of the virtual photon with the u quark, with 
the d quark as a spectator; the bottom diagram represents the 
interaction of the virtual photon with the d quark, with the u 
quark as the spectator. In both diagrams the d is represented 
by a d quark traveling backward in time. 


III. THE TRIANGLE DIAGRAM AND 
INGREDIENTS 


A. The triangle diagram in 
charge-conjugation-invariant GST 


The elastic electromagnetic form factor for a pos¬ 
itively charged 7r+ consisting of a rt and a d quark 
is obtained from the pion current, which—in impulse 
approximation—is the sum of two triangle diagrams, in 
which the photon couples either to the u or the d quark. 
The two diagrams are depicted in Fig. 

In the charge-conjugation invariant formulation of 
CST [5T1 the 7r+ current is given by 


2e 


J{(+(P+, P_) = eF^{Qy{P+ +P_r = i — 


tr 


fkO 


-i - 


tr 


fcO 


r{k,p+)S{p+)j>^{p+,P-)S{p-)T{p_,k)S{k) 

rik,p'_)s{p'_)f{p'_,p'+)s{py)Tip'+,k)sik) 


( 11 ) 


where p± = k + P±, p'^ = k — P±, and j^{p+,p-) is the dressed current for off-shell quarks (defined below). We 
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assume equal masses for the u- and d-quarks, so the u and d propagators are identical. Using this it has been shown 
in our earlier paper |20j that the second contribution to the form factor can be transformed into the first one, and 
the two can be added together. This gives 


'fcO 


J^+{P+,P-) = ie / tr r{k,p+)S{p+)jf^{p+,p_)S{p_)T{p_,k)S{k) 


( 12 ) 


As in Eq. Q, the “fcO” in Eq. 0 stands for the charge- 
conjugation invariant CST prescription of how to per¬ 
form the fep-contour integration. It requires taking all 
quark propagator-pole contributions of the u and the d 
quark into account. The triangle diagram has six propa¬ 
gator poles in the complex fco plane, three positive-energy 
poles in the lower- and three negative-energy poles in the 
upper-half plane. The “fcO” prescription requires averag¬ 
ing over these poles. In the Breit frame, where 

P± = (Po,0,±iQ) 

<Z=(0,O,Q) (13) 


with Pq = + /i the pion mass, and Q the 

photon momentum transfer, the poles of the spectator 
d-quark are locat ed at fco = jyEfe — rjie, with ij = ±1, 
Ek = ^/m? -I- k^, and the poles of the active u-quark 
with momenta p- and p+ are at 


/co = ±E_ - Po T ie , 
ko = ±E+ - Po T ie , 

respectively, where 


E± 




(14) 


(15) 


The active poles can be collectively written as fco = 
— Pq — 77 ie, with ry = ±1 and 77' = ±1. Since the 
square roots in Eq. (15) are positive and p± = k + P±, 
77 = 1 denotes the positive-energy poles, and 77 = —1 
denotes the negative-energy poles of the struck u-quark. 
The index 77' distinguishes between the poles of the active 
quark before and after its interaction with the photon. 


B. Pion vertex function 

One input for the pion form factor calculation is the 
pion vertex function T. Instead of solving the full CST- 
BSE, which will be the subject of another paper, we use 
the approximated pion vertex function 

r(pi,P2) = Goh{pl)h{pl)-f^L{p '^), (16) 

where pi = p + P/2 and p 2 = p — P/2 are the quark 
momenta in terms of the relative momentum p and the 
total pion four-momentum P = (Pq,P'), with Pg = 
\/ P'^; furthermore, h{p1) is a strong quark form 
factor normalized to h(rri^ = 1, (where is the dressed 
quark mass in the chiral limit), and Gq is the inverse 


r 


norm of the wave function. The additional structure 
function L{p^) is a placeholder for the dynamically cal¬ 
culated pion vertex function, to be eventually obtained 
from the solution of the CST wave equation for the pion 
bound state. All we currently know is that L{p^) = 1 in 
the chiral limit, and with this choice the form (16) co¬ 
incides with the one introduced in Ref. [20] and can be 
understood as a finite-pion-momentum extension of the 
chiral-limit pion vertex function in the pion rest frame 


r(p,p) = G^o/*"(p")7' 


(17) 


which is obtained from solving the pion CST-BSE in the 
chiral-limit with a delta-function kernel of the form 

V(p,k) = ^H{pi,p2;k,k2)i2TTfEkS^{p-K){-/f^^-f^) 

- k ) {1^ ® It.) , ( 18 ) 

with C the strength of this interaction (G —>■ in the 

chiral limit m), and with H a shorthand notation for a 
product of four form factors 


H{pi,p2-,h,k2) = h{pl)h{pl)h{m^)h{kl 


(19) 


The first line of Eq. (19) gives the result for PI in the case 


when particle 1 is on-shell in the initial state [so that ki = 
(Pfc, k), kf = m‘^,k 2 = {Ek — Pg, k — P')] and when both 
particles are off-shell in the final state. The second lines 


of the two Eqs. (18) and ( |19| ) give the results in the chiral 
limit in the rest frame, when rug —>■ 0,77i —>■ p —>■ 0. 

Note that a linear confining part of the interaction has 
been omitted from (18); it was shown in Ref. |2T| that 


this part does not contribute to the pseudoscalar bound- 
state equation in the chiral limit. 

From here on, all quantities are assumed to be given 
in the chiral limit, except for the pion mass which is kept 
finite. 


C. Quark current 


In a consistent pion form factor calculation, the dressed 
quark current should also be calculated from solving a 
(inhomogeneous) four-channel CST-BSE. This will, how¬ 
ever, be the subject of a future paper and we use, for 
simplicity, the Ansatze proposed in Ref. |20j . which ap¬ 
plies the framework introduced by Riska and Gross [29] 
to ensure gauge invariance when strong form factors are 
attached to interaction vertices, such as in Eqs. (18) 
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and (19). These form factors can be equivalently moved 
to the quark propagators, leading to so-called damped 
propagators 

S{p) = h^p^)S{p), (20) 

reduced vertex functions 

Tr{pi,p 2) = h~\pl)T{pi,p2)h~\pl) = Goqs , (21) 
and reduced currents 


j^{p',p) = h-\p'^)r{p\p)h-\p^). ( 22 ) 

The reduced current is then required to satisfy the Ward- 
Takahashi identity involving the damped propagators, 

q^jRt.ip',p) = S~^{p) - S~\p '), ( 23 ) 

to preserve gauge invariance. The simplest possible solu¬ 
tion to this equation for pointlike quarks (i.e. quark form 
factors equal to 1 and no magnetic moment) is given by 

[ 201130 ] 

Jr^P^p) = f{p',ph^ + S{p',p)A{-p') 7 '" 

+S{p,p'h^ A(-p) -f g{p',p)K{-p') y'" A(-p). 

( 24 ) 


The functions /, 6, and g, are determined in terms of the 
strong form factor and the mass function through the 
Ward Takahashi identity: 


9iP,p) 

S{p',p) 

fiP,p) 


4MM' - /i'2) 
h'^h''^ (p'^ — p^) ^ 

2M' {M' - M) 

fl/2 ^•p/2 — p2'j ’ 

m 2 - p2 M'2 _ y2 

fl2(^p/2 _ p2'^ / i ' 2 ( p '2 — p'^y 


(25) 

(26) 
(27) 


where we have introduced the short-hand notation h = 
h{py, h' = h{p'^), M = M{py, and M' = M{py. 
The mass function M{p'^), calculated from the chiral 
limit CST-DE using the kernel in Eq. (18), is equal to 
h?{p^), which reduces the off-shell form factors /, d, 
and g to 


9{p',p) 

f{p\p) 


—2S{p',p) —)■ 4to 


2 (fe^ - h'y 

^ (p'2 — p2 ^ 

4<?(P,P)+ f^2i^ay,2_p2) ■ 


(28) 

(29) 


In order to study the influence of the running dressed 
quark mass M(p2) on the pion form factor, we will also 
consider in this paper the case of fixed dressed quark 
masses by setting M{jp) = m^. In this case the off-shell 
form factors reduce to 


9{p',p) 


4m2 (/i2 _ 

(p'2 — p2) 

Hp\p) ^0, 


fip'.p) 


to2 — to 2 — p'2 

h2(p'2 — p2^ h'2(p'2 _ p2'J 


(30) 

(31) 

(32) 


IV. CONTRIBUTIONS FROM THE QUARK 
CURRENT 


Substituting (20), (21), and (22) into the pion cur¬ 
rent (12) gives 


jyip+,py = ie[ G\ 

JkO 


2 h^hlhlL+L_ 


xtr 


kO ^ 

,5 -\j(^ \ Ak' 


1 N{p+)j^{p+,p-)N{py-f A(fc) 


(33) 


where we used the short-hand notation D = M^ — and 
D± = M|. — for the denominators of the propagators 
[with M = M{k^) and M± = M(p^)], and h — h{ky 
and h± = h{py for the strong form factors of spectator 
and active quarks, respectively, and L± = L{py (with 
p± = k + |T±); these additional structure functions are 
assumed to be functions of the relative momentum, but 
they might also depend on p± individually. 


Next we insert (24) for the the quark current into (33), 


which gives three contributions to the pion current asso¬ 
ciated with the off-shell form factors /, 6, and g: 


Ji = 4’^ + 4^ + . 


(34) 


In the following we analyze these three contributions 
separately and calculate for each of them the contribu¬ 
tions from the spectator and the active poles. 

Note that, from the off-shell structure of the g{p',p) 
term in Eq. (24), we can already anticipate that only the 
spectator pole will contribute in this case. This illus¬ 
trates one of the advantages of breaking the calculation 
of the form factor into these three terms, instead of only 
concentrating on the total contribution of each of the six 
poles. The decomposition given by Eq. (34) turned out 


to be very useful for cross-checking our analytical and 
numerical calculations. 


A. /-term contribution 

The / contribution to the pion current is 

g: 


where 


Q2^+A^-A+L-f{p+,P-) 


/fco 


DD^D 


xtr 


= le 


iV(p+)7'^V(p_)A(-fc) 

^h^L+L^A{k+,k-)Aff 


G^ 


'kO 


{pl-pl)D 


A(p+,p-) = 


hi 


Ml -pI M^-pi' 


(35) 


and the trace is 


(36) 
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J\ff = 4fc^(p_ • p+ — M_M+) + (MM+ — k ■ p+) + Ap^lMM- — k ■ p-) 

= Ak>^{MM+ - M_M+ + MM_ - k ■ p-+ p-■ p+- k ■ p+) + AP^{MM+ + MM_ - k ■ p_ - k ■ p+) 
+2q^{MM_ - MM+ + k ■ p+ - k ■ p_). (37) 


Note that p\ = (/co + -Fb)^ — (k± ^q)^, so that p\ —p^_ = —2kzQ, and excluding the factor the integrand is even 
in k. The term in Af-^ is odd in kz and hence this term integrates to zero - a consequence of current conservation. 
For the k^ term in we introduce 2Pq = {P+ + P-)^ and use the relation k^ —>■ P/^ko/Pg, which holds since the 
rest of the integrand is even in k. Then, the remaining terms reduce to 


Aff = AP^^{MM+ - M_M+ + MM_ - k^ + P+■ P_) + AP^{MM+ + MM_ - 2k^ - 2k ■ Pq) 


= AP,^ 


(MM+ + MM_) 1 + — - 1 + 


2Pn 


+ (P+ • P_ - M+M_)^ - 2koPo 


= 8P^ <! M 


' M+ + M_ 


- + [M{M+ + M_) - M+M_ - fc^l — 

2Pq 


fj^ 

2Pn 


(38) 


The / contribution to the form factor then becomes 


Fim = -iGi [ 

J kO 


A(p+,p_)x(fco, Q) 
kzQPo 


,(39) 


where 


x{ko, Q) = ^ {Po [M{M+ + M_) - 2k^] 

+ko [M{M+ + M_) - M+M_ -k^ - p^]} . (40) 


For fixed dressed quark masses, M = M+ = M_ —>• 
and x{ko,Q) reduces to 


x(fco) Q) 2Pq + fco 




(41) 


Note that in the chiral limit p = 0, such that the contri¬ 
bution of the spectator in the third term of this equation 
vanishes, and, in this limit, the entire result comes from 
the active quark poles in A. Next we perform the kp 
integration. 


1. Spectator pole contribution 


The relative momenta that appear as arguments of the 
L functions become 


Pu'v — (^ + 2^n') — k"^ + k ■ Pr,' + ^ 
= '^x + ^ + pEkPp - p'kz j ■ 


(44) 


In (421 we have also used the abbreviations 


hrj'ij = , 

Mjjijj = , 

Lri'n = L{Pri'n) j 


(45) 


and 


Xr)(Q)= x{vEk,Q) 

- 2m2] 

+ 7;V [ntxi^+,V + ^-,v) - - ^x ~ ' 

(46) 


First we calculate the spectator pole contribution to 
F^m in ( |39[ ). Averaging over the poles at ko — rjEk 
gives 


Fi’%Q^) = Gl 



P XriiQ)krfinF+,riL— 

2PpkzQ{M^,^-p^z,'u) 


(42) 


where p^,^ is the square of the active quark momentum, 
p^,, evaluated at the spectator pole at kp = pEk^ such 
that 

Pu'v = + ^o)^ “ (^ + 

= rn^ + + 2pEkPp - p'kzQ ■ (43) 


For fixed quark masses, Xn(Q) ~PP ^and the con¬ 
tribution to the form factor simplifies to 


F/’^(Q2) ^ -Gl 


d^k 


P 


(27r)3 APpkzQ 

V- PP'K' 


2 r r 


7 ], 7 ]' 


K-Pv'v 


(47) 


In Appendices |A 2 c| and |A 2 d| we analyze, for fixed quark 
masses and L{p‘^) = 1, the small and large behavior of 
the spectator contribution to the / term. We find that 
the contribution has a Q~^ falloff at large consis¬ 
tent with the experimentally observed flattening of the 
Q'^F^{Q‘^) curve. 
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2. Active pole contributions 


Next we calculate the active pole contributions by aver¬ 
aging over poles in the upper half plane at kg = —E± — Pq 
and in the lower half plane at ko = E± — Pq . These cases 
are described by the subscripts {rj'r]} in the general ex¬ 
pression ko = rjErf' — Pq- We need to evaluate, at the 
active particle pole {rj'r]}, the squares of the spectator 
momentum, fc, of the active quark momentum, Pri", and 
of the relative momentum, p^//. They are denoted, re¬ 
spectively, by 


kl^ = {rjE^, -Pof- k' 

= ml + - 2riEr,'Po + p'kzQ + 

^ 1 

(Pri'^v'n ~ 

iP-r,')r,'ri = {vEr^'f “ O) 

= m\ + 2r]'k:,Q = 

iPu')ri'v = {v^v’ - \Po) -kl- (k^ + \v'Q) 

^72 . ^ 2 ^2 

- 2^v'v + 2^x ~ 

iP-v')v'n — {yEri’ — 2^^'} ~ ^-L ~ ~ 4^ 


— iPv')v'v + p'kzQ ■ (48) 

Since {p'^f)r)'r] = we need only a simplified notation 
for namely 


M,, (49) 

L+L_ always occurs in the form of a symmetric product, 
such that 


Efi'r] — E[{Prji)n'v\k^[iP--n')n'n\ ■ (^ 9 ) 


We obtained the somewhat unexpected result that the 
/-term contribution from the active quark poles has the 
same Q~^ asymptotic falloff as the one from the specta¬ 
tor quark pole. This finding is independent of the precise 
functional form of the quark form factors and of details 
of the pion wave function. 


B. (5-term contribution 


Next we look at the contributions from the S terms in 
the quark current. Because the S form factors are propor¬ 
tional to the difference of the mass functions, there are 
no contributions to this term in the case of fixed quark 
masses. For running quark masses the (5-term contribu¬ 
tion is given by 


= i 


le 


Glhlhlh^L^L_ 


tr|iV(p+) 


Jko DD+D- 

s{p+,p-)H-p+h^ + s{p-,p+h^H-p-) 


xN{p_)N{-k)'^ 
'''J.o 2kzQD ^ 


hlAfl hlAfi ] 


(53) 


where 6M = M+ — M_, and the traces are 


= tr{7'^(M_ -t ^_)(M - 1^)} 

= 4/c^(M - M_) -k AP^M + , 

A7^=tr{(M++|i^)7'^(M-|^)} 

= Ak>^{M-M+)+AP^M-2q>^M. (54) 


As before, the coefficient of the term is odd in kz, and 
thus integrates to zero, a consequence of current conser¬ 
vation. Since the factor multiplying k>^ is even in kz, we 
can substitute k^ —)■ koPg /Pq in the trace terms, which 
gives 


At the pole indicated by {p'p}, x{ko, Q) is 


Xn'viQ) = xivEr,' - Po,Q) 




M^(^(M^( -k m^) - 

"'v'v 

+{vE,^, - Po)[M,j'^{Mn' + m^) - - k^,^ - E 


!}■ 

(51) 

The /-term contribution of the active poles to the form 
factor then becomes 

Fpim = 011 . (52) 


v,v' 


AkzQEzi' Pq 


A2a 


and 


A2b 


we 


where hrj>rj = h{k‘^,^). In Appendices 
analyze the small and large behavior of this contri 
bution for fixed dressed quark masses and L(p^) = 1. 


Mi 


Mt 


APi 


APi 


{M - M+)— + M 
Eo 

{M - M_)^ + M 
Eo 


Fim = -i f 

Jh 


GIML+L 6M 


kO kzQPo M2 - *2 
X {Ei{p){Eo + ko) - j: 2 {p)ko} 


where 


Ei{p)=M 


hi 


Ml-pi Mi- pI 


, , M_M M+hl 
S2(P)= .., ~, + + + 


Ml-pi Ml- pI' 


(55) 


Then the contribution to the form factor becomes 


(56) 


(57) 


It is explicit in the three last equations that both active 
and spectator quark poles contribute to the 6 term of the 

























pion form factor. Next we perform the ko integration. 
Averaging over the spectator poles at ko = rjEk gives the 


spectator contribution 


= Glj 




(27r)3 Ak.QPoEk 






M+^n - M_^n 
M"^, — p^, 

n'r) t'n'r) 


[m^Po + pEk{m^ - Mrj'n)] , (58) 


and averaging over active quark poles at ko = pE^ji — Po gives 

d3fc 1 


= -Gl 


(27r)3 Ak.QPo 


V ^r]'rj^7]'rj ( Mjji — TTl^ 


A, 


v' 




[m^Po + pEjj> 


(59) 


r 


C. gr-term contribution 

The remaining contribution comes from the fully off- 
shell g term in the quark current, which reads 


strong quark form factors of the simple form 

2 

/ A ^ m ^ \ 

Kp^) = 


^2 — ™2 




(63) 


J£’^ = ie 


G^h^hihiL+L. 

-^- -9 [P+,P-) 


IkO 


DD+D 


xtr|lV(p+)A(-p+) 7 ^A(-p_)iV(p_)A^(-fc)| 
2Glh‘^{h\ - h^_)L+L_ki^ 


= —le 


' kO 


k,QD 


(60) 


There is no contribution from the active quark poles to 
this term, because the denominators and D_ can¬ 
cel with N{p+)A{—p+) and A(—p_)A^(p_), respectively. 
Hence, the only contribution comes from the spectator 
poles. Excluding the factor of k^, the integrand is even 
in k. Hence, we can replace k^ koPo /Po and the 
( 7 -term contribution to the pion form factor becomes 


Fa{Q^) = -iGlf ■ 

J kO 


h^{hl - h'i)L+L_ko 
k^QDPo 

Then, averaging over the poles at ko = pEk gives 


(61) 


F^^m =Gl I 


d^k 1 
(27r)3 dk^QPo 


'y ^ VP ^ri'ri^+,V^-,V ' 


ri,ri' 


(62) 

It is interesting to observe that the ( 7 -term contribution is 
identical for fixed and running quark masses, because the 
depe ndence on the m ass function cancels out. In Appen¬ 
dices Ala and A 1 b we analyze, for L{p^) = 1 , the small 
and large behavior of this contribution. Similarly to 
what we found for the /-term, the ( 7 -term contribution 
also falls off as and is also independent of the func¬ 
tional form of the quark form factors. 


V. RESULTS 

In our previous paper. Ref. m, on the pion form factor 
in the relativistic impulse approximation, we have used 


Here A^ is a cut-off parameter. 


Both, and A^ are 
determined by a fit of the quark mass function M{p^) = 
m^h 


^(p^) at negative p^ to the lattice QCD data [3T] 
extrapolated to the chiral limit m- The fit gives Ay = 


2.04 GeV and my = 0.308 GeV. Note that h(p^) of (|63|) 


has a pole at A^ = P^- Our previous calculation only 
explored the values of p\ at the RIA spectator particle 
pole (at ko = —Ek), which, from Eq. (43), is bounded by 


p2, _ = .^2 ^2 _ 2EkPo - p'k^Q 

< -I- < A^ . 


(64) 


The contributions from the active particle poles and the 
spectator pole at fcg = E^ included in the present C-GIA 
calculation will probe the structure of the form factor h 
also at large positive p^ and will therefore depend on the 
definition of h in this region. To study this sensitivity, 
we adopt a generalization of the piecewise form proposed 




h{p^) = 


Al -p 2 


AA(a) 


Q!^A^ — 

(a^A^ + p'^ — 2 s+ 



where s+ < A^ is some fixed value (given below), and 


the normalization factor 


AA(a) = 


(A^ - m^){a^A‘^ - s+) 
(q;2A| — to|)(A 2 — s+) 


n 2 


( 66 ) 


makes h{p‘^) continuous at = s+. Note that this def¬ 
inition of h(p^) in the region p^ > ensures that it 


has no pole at p^ = A^, or anywhere in the region if 
a^A^ > s+. Varying the value of the parameter a allows 
us to study the sensitivity of the pion form factor to the 
functional form of h(p^) in this region (where it is not 
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FIG. 2. (Color online) The chiral limit mass function M{p^), 
for a = 3 (dashed, brown), q = 1 (solid, orange), and a = 0.5 
(dotted, purple) compared with the lattice data [ST] extrapo¬ 
lated to the chiral limit. 


constrained by lattice data). If a = 1, the new defini¬ 
tion (65) is identical to the one proposed in [2T] , and we 
emphasize that when < s+ , the form (651 is identical 
to Eq. (63) used in our previous calculations of the pion 
form factor in the RIA [5T] and, therefore, the results of 
these calculations remain unchanged by using ( [^ . We 
choose s_|_ = (A^ — rriy^Y/A = 0.752GeV^. For a = 1, 
hljp') is symmetric about = s+. Figure shows the 
mass function M{p^) = m^h^{p^) for a = 0.5,1, and 3 
together with the lattice data extrapolated to the chiral 
limit. 


We have calculated the form factor for both, fixed 
dressed quark masses by setting M{p^) —)■ and run¬ 

ning dressed quark masses, with M{p^) = m^h?'{p^). In 
both cases we keep the strong quark form factors in the 
pion vertex functions and in the quark current as defined 
by Eq. ( [^ . All results presented here were obtained 
with L{p^= 1. 

In Fig. we compare the spectator and active 
quark pole contributions, calculated with running quark 
masses, for different values of the pion mass, namely 
pL = 0.14, 0.42, and 0.6 GeV. Also shown in this fig¬ 
ure is the asymptotic function \Q~^ which gives a good 
fit to the pion form factor at large as can be seen in 
Fig.i Gomparing the slopes of this asymptotic form to 
the form factor curves shows that all of the contributions 
have the correct 1/Q^ falloff at large . One notices 
in the results of Fig. that, as the pion mass decreases, 
this asymptotic behavior sets in for smaller values of . 
Note also that both the spectator and active pole contri¬ 
bution are normalized by the same factor to ensure that 

f,^(0) + f;(0) = i. 

Figurej^shows the ratio for fixed and 

running quark masses, and different values of p. It is a 
surprising result that these hxed and running mass cal¬ 
culations are very close to each other over the entire 
range, with the possible exception for large p and large 
where the running mass results are larger by about 
10 %. 


Cf (GeV^) 



(GeV^) 



FIG. 3. (Color online) Comparison of the normalized specta¬ 
tor pole (red, dotted) and active pole (blue, dashed) contribu¬ 
tions, calculated with running quark masses, for ^ = 0.14 GeV 
(top panel), for p = 0.42 GeV (middle panel), and fj, = 0.60 
GeV (bottom panel). In each panel the sum of active and 
spectator pole contributions is the solid purple line which 
approaches 1 at small . The asymptotic function 
which fits the high form factor data, is shown for compar¬ 
ison. 


We have also looked at the sensitivity of the pion form 
factor to the shape of the h form factor in the timelike 
region, where it is not constrained by lattice data. 
In particular, we have varied the parameter a between 
0.5 and 3. We find that the pion form factor is quite 
insensitive to these variations; only the active quark pole 
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FIG. 4. (Color online) The data |32H40 | for the pion form 
factor shown against the simple model F^{Q'^) = (solid 

line). 



FIG. 5. (Color online) The ratio F^jF^ for fixed (dashed 
lines) and rnnning (solid lines) qnark masses, and different 
values of fi. The pairs of curves, from top to bottom are the 
results obtained with fi — 0.60 (brown), 0.42 (orange), and 
0.14 GeV (purple). 


contributions at large display any sensitivity at all. 
Figure]^ shows the results for ^ = 0.42 GeV. 


VI. SUMMARY AND CONCLUSIONS 

The present study of the pion form factor using the 
Covariant Spectator Theory (GST) improves on our pre¬ 
vious work [20) . which was done in the relativistic im¬ 
pulse approximation (RIA), where only one quark pole 
contribution to the triangle diagram for the pion cur¬ 
rent, the negative-energy spectator pole, was calculated. 
In the present work we take all six quark pole contribu¬ 
tions into account. These are the positive and negative- 
energy spectator poles and the positive and negative- 
energy poles of the quark that interacts with the photon. 
The latter are referred to as the active poles. The inclu¬ 



Cf (GeV^) 



FIG. 6. (Color online) The pion form factor contributions 
calculated for a = 3 (dashed, brown), a = 1 (solid, orange), 
and a = 0.5 (dotted, purple), with p = 0.42 GeV and with 
running quark masses. The top panel shows the spectator 
and the bottom panel shows the active pole contributions. 


sion of both, positive and negative-energy poles is nec¬ 
essary in order to preserve charge-conjugation invariance 
in GST [^121], hence we refer to the present calculation 
as the C-symmetric complete impulse approximation (G- 
GIA). 

For the calculation of the pion form factor in C-CIA we 
use the same off-shell quark current as introduced in [20] . 
It satisfies the Ward-Takahashi identity and therefore our 
pion current is conserved. Our results for the pion form 
factor in C-CIA depend also on the pion mass fi, similar 
to what was found in our previous RIA calculations [5D] . 
As expected from the pole analysis therein, we find that 
for small ^ the active quark contributions are as impor¬ 
tant as the spectator contributions, over the whole range 
of Q^. For large /x and large Q^, the active pole contri¬ 
butions are smaller than the spectator contributions by 
about 30%. For small /i, the spectator and active pole 
contributions are nearly identical, not only in magnitude 
but also in shape, even for large Q^. This is a surprising 
new result. Nevertheless, it confirms the correct large-Q^ 
behavior obtained in |20j . which remains unchanged by 
the inclusion of the active poles, and sustains the flatten¬ 
ing of Q^F)r((5^) for large Q^. 
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We emphasize that the pion vertex function, as intro¬ 
duced in Ref. [20], is a very simple Ansatz, that can be 
understood as a finite-pion-momentum extension of the 
pion vertex function in the chiral limit. The similarity 
between spectator and active pole contributions may be 
the result of using this simple pion vertex function with 
L{p^) = 1. 

We have also studied the effect of the running dressed 
quark masses on the pion form factor. We find that for 
small p the results for the pion form factors calculated 
with running and with fixed dressed quark masses are 
nearly identical, and only for larger p. and we observe 
some differences, which are, however, quite small. 

Related to this we have also investigated the sensitivity 
of the pion form factor on the strong quark form factors. 
In the C-CIA, the quark form factors, and accordingly 
the quark mass functions, are tested over the whole 
range. In particular, at positive where their shapes are 
not constrained by lattice QCD data, this offers the pos¬ 
sibility to obtain important information about the mass 
function in this region. However, we find that for small 
pion mass the shape of the pion form factor, even in the 
small region is quite insensitive to the functional form 
of the strong quark form factors. The early onset of the 
Q~^ behavior for small p may be connected to this ob¬ 
served insensitivity. Only the active poles for large p and 
show a small dependence. 

It should be mentioned that in this first C-CIA study 
of the pion form factor, we were rather constrained in 
our choices for the pion vertex function. In particular, at 
small pion masses our form factor persistently falls below 
the data at small . A more stringent comparison of our 
calculations with the data will be possible once precise 
pion vertex functions as solutions of the exact CST-BSE, 
as well as a dynamically calculated fully-dressed quark 
current, become available. Work on both subjects is cur¬ 
rently under way. 
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Appendix A: Small and large Q behavior of the pion 
form factor 


1. gr-term contribution 


Here we investigate the small and large Q behavior of 
the g-term contribution Eq. (611. 


a. Small Q behavior 

It is easy to obtain the small Q behavior if p ^ 0 (the 
form factor diverges for all Q when p = 0). Then, to 
order Q, 

dh'^ 

(Al) 

where 


^ , (A2) 


and 


S = HPv) ■ (A3) 

(Here symbols with “hats” should not be confused with 
on-shell four-vectors.) Note that —)■ poo as fc —>• oo, 

only if p 7 ^ 0. Since does not depend on ry', the p' sum 
contributes a factor of 2, and the small Q limit becomes 


^im F®(Q2) = -Gl I 


d^k 1 dh'^ 

(27r)3 2p^^dp‘^' 


(A4) 


where we have used that Pq p as Q —>■ 0. 


b. Large Q behavior 


At large Q the arguments of h can only be small (and 
hence the integrals large) if k^ is very large (so that k^ ~ 
Ek). Expanding the arguments for large kz and Q gives 


Pv’v ^rn^ + p^ + Q 

2pp^\k, 




2\k, 


E- 


Q 


= pIc, ■ 



(AS) 


For to be small at large kz requires that the sign of 
kz be chosen so that p\kz \ — p'kz = 0, which will be true 
only for positive if ?7 = p', and negative kz ii p = —p'- 
Hence the integrals for p = p' run from 0 < kz < oo, 
and for p = —p' from —oo < kz < 0, and with this 
understanding, the arguments go to 


-I- -I- pQ[ 

m‘^+ p'^ - pQ[ 


2fc^ 


2fi^kz 1 

’ 

1 

> 


p = p' 
p = -p'. 

(A6) 
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We obtain 


Rs 

lim F^Q^) = — (A7) 

3->oo 


where, transforming kz = Qk'z, 



where we introduced the convenient notation 

and, in the second equation, we changed the sign of k"^ 
allowing the two terms to be combined into one (with a 
factor of 2). Note that the singularity that was at kz = 0 
is now suppressed by the property h(±oo) = 0 , and that, 
because fi y^O the integrals converge at both k” = 0 and 
k” = oo. The integrals peak at 


+ k 


peak 


2/i 


(AlO) 


Notice that if 77 = —1, and k” = 1/2, the argument 
becomes 


-Foo,— 


= -k 


± ■ 


(All) 


2. /-term contribution 


Here we investigate, for fixed quark masses, the small 
and large Q behavior of the /-term contribution F/ (Q^), 
Eq. 


a. Small Q behavior of active pole eontributions 


The small Q behavior can be obtained by expanding 
near Q = 0 to order Q, 


1,2 12 , 

hv'v-h-v+ Ek 


'dp/ 


(A12) 


which is nonsingular as long as 2m^ > p, and taking the 
y —)• 0 limit in the rest of (52) gives 


^imF/’“(Q2)=G2 J 


d^k 1 

(27r)3 p 

viyEk - /i) dh'i^ 

^ 2r]Ek - fJ. dp/ 


(A14) 


where the two identical t]' terms have been combined, 
giving a factor of 2. Changing ry — > —77 in the sum gives 
a result which reduces to Eq. (A4) if p = 0, 


^im F/’“(Q2) = -Gl J 


d^k 1 

{2tt)^ P 


sivEk + p) dh^ 

^ 2-qEk -f p dp 2 ■ 


(A15) 


b. Large Q behavior of the active pole contributions 


The large Q behavior of (521 follows from arguments 
similar to those used for the study of F®. Anticipating 
that the dominant contributions come from large kz, we 
expand the arguments as follows: 


1 


1 


kz, ~ + p^ - r]Q\kz + xJf'QI + 'n'Q kz + -p'Q 


-vQ 

= kt.. 


mt 


ki 


2|fc. 


\v'Q\ 


2p^\kz, + ip'QI 


(AI 6 ) 


Since Q can be chosen to be positive, the argument is 
finite at large Q only if 


'n\kz + \^Q\ -v' kz + ^v'Q 


= 0 , (A17) 


which leads to the conditions 


kz>-^vQ v = v' 


(A18) 


and 


kz < v = -v' ■ 


(A19) 


Because of these conditions, when 77 = 77 ' it is convenient 
to shift the kz integration to kz ^ k'^ — \riQ, so that 
kz + s-iid the region of integration is k'^ > 0, 

while if 77 = —rj' the shift is kz —>■ k'^ + \riQ, so that 
kz — ^vQ and the region of integration is < 0. 

With these transformations, the arguments become 


where p^ and were defined in Eq. (A2|. Using 


= -JLL. , (A13) 

2r]Ek - p 2riEk - fJ. 


^00 


777^ + p^ - r]Q[ 
777^ + p^ + r]Q[ 


2k'^ 


2k' 


2/x^A:^ 1 
2fi^k'^ 1 


v = v' 
v = -v' ■ 
(A20) 


Xv'vi^) = vEk + fk 
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Note that 

fcL =ploiv^ -V, P -V, 


K) ■ (A21) 


This means that the results obtained from the study of 
the (;-term can be used here: the high form factor 
is dominated by > 0 when rj = r]' and < 0 when 
r] = —rj' . In these regions, after shifting kz —> k'^ in the 
original definition (51) (as discussed above), the large Q 


behavior of Xr)'ri becomes 

X^'viQ) - vlK 


1 2fi^Q\k'J{r]\k'J - kQ) 

2^ Q^V{7j) 


Q^V{r,) L 


Q^{ml + kl)m 


1 


vQ) 




+A^^X\\K\ - Iv) 
= Qxr{v, \K\) > 


(A22) 


where, in the second from last line, we anticipate the 
substitution k'^ = Qk'^ to be made below, and 

Q^V{r]) = vQ^irnl + kl) + A-qix^kf - 2^^Q\k',\ 

= ri{m\ + k\) + Aruj^k'^"^ - 2^^|/c"| (.A23) 

Note that the reduced xr{v^ l^il) = XRiv^ k") if k'l > 0, 
and xk(? 7 , |fc"|) = —fc") if fc" < 0. Furthermore, it 

has the symmetry XRiv, k") = Xr{-V, -k”)- 


Finally, at large Q the result (52) becomes 


where, making the substitutions discussed above, 
d^A;L 


(A24) 


= Gl 


(27r)' 


■E 


V 


/■“ XXRiV,k'^) 2 r //2 1 

Jo 2k'l{k'^-\v) 


= Gl 


/-oo 2|fc"|(fcg 

^2 




[ f°° XxRiv,k'') ^2, //2 1 

J (2XX k'X'l-h) ^ 

(A25) 


where, in the second line, we have combined the two in¬ 
tegrals over k" by changing k" —>• — fc" in the second 
integral. When 77 = 1, the integral has a principal value 
singularity at fc" = f, which can be easily integrated 
over. 

c. Small Q behavior of the spectator pole contribution 

Evaluation of this contribution follows the discussion 
of the g-term contribution. We obtain immediately 


lim = a 

Q-s-O ^ ' 


= Gn 


d^k fi 
(27r)3 2 

d^k 1 


E 


d r 


dp2 Lto 2 _ p2 

E " 


(27r)3 2fi ^ 2r]Ek + H 


K 


dhl 


P- 


'' f. (A26) 


y{2'qEk + li) dp. ^ 

The value of the total form factor at Q = 0 is obtained 
by adding this to (A4| and (A15), giving 


^imF^(Q2) = GlJ 


d^k 1 

(27r)3 2/r 




- 2 - 


E^ 

77 

dhl 


I 

I {2pEk + X dp: 


(A27) 


d. Large Q behavior of the spectator pole contribution 

The large Q behavior also follows from the discussion 
of the p-term contribution. Using (A 6 ) and repeating the 
steps leading to (A 8 ) gives 




(A28) 


where, transforming k^ = Qk", and using the definition 

d^, 


d/,s ^ _pi2 f d^fc-L \ " j dfc" 

Oj (2^)3 2 Jo k'f ml-p'X 

U! ?,2r„//2 


/■° dfc" hXXn\ 

/_oo k'l 


= 2Gi 


d^fci f 

xr Jo 


dfc"E 




00 , 7 ]] 


rrif -|- fc^ -I- 2p?k'f ( 2 fc" -|- 77 )' 


(A29) 


Combining the results for and i?® gives the 

total coefficient for the Q~^ falloff. 
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